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It is investigated theoretically that magnetized Bose-Einstein condensation (BEC) with the 
internal (spin) degrees of freedom exhibits a rich variety of phase transitions, depending on the 
sign of the interaction in the spin channel. In the antiferromagnetic interaction case there exist 
always double BEC transitions from single component BEC to multiple component BEC. In the 
ferromagnetic case BEC becomes always unstable at a lower temperature, leading to a phase 
separation. The detailed phase diagram for the temperature vs the polarization, the spatial spin 
structure, the distribution of non-condensates and the excitation spectrum are examined for the 
harmonically trapped systems. 

KEYWORDS: Bose-Einstein condensation, Non-condensate, Gross-Pitaevskii equation, Bogoliubov equation, Crit- 
ical Temperature 
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§1. Introduction 

There has been much attention focused on the Bose- 
Einstein condensation of alkali-metal atom gases both 
experimentally and theoretically recent years since the 
experimental realizations in 1995. 1_3 ) Interest ranges 
widely from atomic physics, quantum optics to many- 
body physics. One of the focused points from the many- 
body physics viewpoint lies in the fact that the Bose- 
Einstein condensation with internal degrees of freedom 
can be realized, thus allowing us to investigate the inter- 
play between superfluidity and internal degrees of free- 
dom. Namely, in the recent experiments on Na atoms 
via an optical trapping method the three hyperfine states 
with mp — 1,0, —1 of F = 1 for the case of Na simul- 
taneously are Bose condensed as named a spinor BEC. 4 - 1 
This contrasts with much widely used magnetic trapping 
experiments where the freedom due to various hyperfine 
states is frozen and BEC is described by a scalar order 
parameter. 

Prior to the realization 5 - 1 of a spinor BEC, only a few 
examples of the superfluid with internal degrees of free- 
dom such as spin or orbital degrees of freedom are known 
like in a neutral Fermion system: 3 Hc 6 ) or a charged 
Fermion system: An advantage for consider- 

ing dilute alkali-metal atom gases lies in the fact that 
these are weakly interacting Bose systems, thus mean 
field treatment of Bogoliubov theory is well applicable, 
which is a controllable approximation. This feature is 
absent in 3 He and UPta. This allows us to construct a 
theory from a microscopic Hamiltonian without any ad- 
justable parameter. This is one of unique properties in 
the present alkali-metal BEC systems. 



In this paper, we consider the phase diagram of the 
BEC with internal degrees of freedom characterized by 
having the atomic hyperfine state with F = 1. Thus the 
triple degenerate components (mp = 1,0,-1) are avail- 
able and can simultaneously Bose condense by an optical 
trapping. It is well known that the condensation temper- 
ature decreases when the spin degrees of freedom exist 
than when the spin degree of freedom does not exist, 
namely the degeneracy lowers the BEC critical tempera- 
ture. Therefore a magnetically trapped system has gen- 
erally a higher BEC transition temperature than an op- 
tically trapped system, if the other conditions are same. 

On the other hand, the total magnetization of the 
system is preserved in BEC because the system is iso- 
lated and the magnetic dipole interaction is weak. In 
this case with mp = 1,0,-1, the two phase transitions 
are expected rather than three when the system has a 
finite magnetization. The majority spin component, say, 
mp = 1 whose density is larger than the others first 
Bose condenses upon lowering temperature T. The sec- 
ond transition occur only when all three chemical poten- 
tials, which were spaced equally by the Zeeman effect in 
the normal state, coincide. Thus the third transition is 
never realized for a finite magnetization from the general 
ground. This double transition phenomenon is similar 
to the split transitions to the Al phase and A2 phase 
of superfluid 3 He-A under the magnetic field. 6 ) It is an 
interesting problem how these transition temperatures 
actually varies when taking into account the polariza- 
tion as a controlling parameter of the system. It is also 
interesting to see how the interaction makes the nature 
of the phase transition alter and the critical tempera- 
tures change relative to ideal Bose systems. The spin 
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dependent interaction constant which characterizes the 
spinor BEC could be cither ferromagnetic and antifcr- 
romangentic. The former is anticipated in 87 Rb 4 ) while 
the latter is realized in 23 Na with F = 1. We will see 
that the sign of this interaction channel is crucial in de- 
termining the low temperature properties even though 
its strength may be small in actual systems. 

In the next section, the phase diagram in a system 
of non-interacting Bosons is discussed as a preliminary. 
A microscopic formulation based on a mean-field the- 
ory of Bogoliubov approximation is introduced and the 
Popov equations which describe interacting Bosons with 
the spin degrees of freedom are derived in Sec. 3, fol- 
lowing our series of papers 8, 9 ) and others 10 -* on spinor 
BEC. Using them, we discuss the phase diagram of the 
interacting Boson systems. The effects of the interaction 
both for ferromagnetic and antifcrromagnetic cases are 
discussed there. The last section is devoted to summary 
and discussions. 

§2. Ideal Bose systems 

In order to facilitate the later discussions on the inter- 
action effects of the spinor BEC, we first treat a system 
of non-interacting particles with F = 1. The system is 
assumed to have a given polarization M and cylindrical 
symmetry along the z direction. It is trapped harmon- 
ically with the trap frequency u)/2tt in the radial direc- 
tion. The three kinds of the particle number N of non- 
condensed atoms, which are specified by the hypcrfinc, 
or spin states vnp = 1, 0, —1, are written as 



N(P\m\) 



E 

n x ,n y =0 



dk z 



exp0(e(n x ,riy,k z ) + \m 



1' 
(1) 

where the chemical potential Hi{< 0) corresponds to each 
spin state (i — 1, 0, —1) and they are split equally by the 
Zeeman effect: /ii + = 2/io to maintain the given 
polarization. The eigenvalue is given by 
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= 1/ksT and L is the length along the z direction. 
N(0Hi) is simplified to 



C = y^k 
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(3) 
(4) 

(5) 



Let us consider the critical temperature T cl , which is 
determined by the condition fi-^ = where the conden- 
sate with the spin state i = 1 becomes non-vanishing 
from high temperature. Since jii = 2/io, the parti- 

cle numbers of the spin state of +1, and —1 are given by 
N(0),N(a) and N(2a) respectively where the dimension- 



less implicit parameter a is proportional to the Zeeman 
splitting. The particle number N and the polarization M 
which are given parameters of the problem are expressed 
in terms of the implicit parameter a as 

N = CT*{ 2 {f(0) + f(a) + f(2a)}, (6) 

M = CT 5 J 2 {f(0)-f(2a)}. (7) 

When all particles in the system are in the state i = +1 
and the fully polarized case: M/N — 1, we denote the 
critical temperature as Tq. This corresponds to the limit 
where only the single component condenses, that is, a — > 
oo. The total particle number at T — T is 



N = CT 5/2 /(0). 
Comparing N in eqs. (6) and (8), we obtain 



Tci 
T 



/(0) 



2/5 



(8) 



(9) 



_/(0) + /(<») + /(2a) 

The relative polarization M/N is derived from eqs. (6) 
and (7) as 



M 
TV 



/(0) - /(2a) 



/(0) + /(a) + /(2a)- (10) 

Equations (9) and (10) constitute the coupled equations 
to determine T c i/T$ as a function of M/N through the 
parameter a. 

Let us consider the second critical temperature T c2 
at which all three components simultaneously condense, 
namely, at which a = 0. The particle numbers of 
the spin components with +1,0 and —1 are given by 
N c + N(0), N(0) and JV(0) respectively, where N c is the 
number of atoms in the condensate. The total atom num- 
ber and the polarization are expressed as 



N 
M 



N c + 3CT^ 2 f(Q), 



N, 



thus giving 



N-M = 3CT c 5 2 /2 /(0) 



(11) 
(12) 

(13) 



We can consider the point which yields same polarization 
M on T c i line as a function of M. The value N - M is 
evaluated from eqs. (6) and (7) as 

N-M = CT*{ 2 {f(a) + 2f(2a)}. 



Comparing eqs. (13) and (14), we find 

T r .2 //(a) + 2/(2a)^ 2/5 



3/(0) 



(14) 



(15) 



Equations (9) and (15) determine the T c i/Tq and 
Tci/Tq as a function of M/N. It is noted that at 
M/N = 0, T c i/T = 3~ 2 / 5 = 0.644. This is reason- 
able because the critical temperature decreases as the 
degeneracy of the system increases. Note that the above 
factor 3 comes from the triple degeneracy. 
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Figure 1 shows the two critical temperatures T c \ and 
T C 2 as a function of the given polarization. It is found 
that at a fixed polarization M/N(^ 0) upon decreasing 
T the system always shows the double transitions from 
the single spin component condensation; the A phase at 
T c i to the three component condensation; the B phase 
at T C 2- This double transition feature is similar to the 
Al and A2 transitions in the A phase of the superfluid 
3 He under an external field. 




-.- M/N 



Fig. 1. The phase diagram of the non-interacting system. The 
axes are the normalized polarization M/N and the normalized 
temperature T/Tq. The A (B) phase is characterized by the 
single (three) spin component BEC. N denotes the normal state. 



0,-1) and expresses the magnetic field acting on the spins 
as the Zeeman effect. The harmonic confining potential 
is V(r) = \mJ 2 r 2 . 

Following our previous paper, 8 ) and the usual proce- 
dure for BEC system without spin, the Gross-Pitaevskii 
equation including the non-condensate as the mean field 
is obtained as 

[h{r)5 i3 - Bij 

+9n!^{\<t>k\ 2 + Pkk)Sij+P* ij \ 
a kl 

+ {Fa)ikP%l(Fa)lj}]<l>j = 0. (20) 

The non-condensate density 



Pki = {4>k (r)ipi{r)) 



(21) 



is introduced. 

The non-condensate density is obtained through the 
excitation spectrum, which is the solution of the Popov 
equations given by 

Y^{AijU q {r,j) - B ijVq (r,j)} = e q u q {v,i), (22) 

3 

Y,{Bfj u q (r,j) - A* jVq (r,j)} = e q v g (v,i). (23) 



§3. Formulation for interacting Bosons 

In order to treat an interaction Boson system, we give 
the formulation based on Bogoliubov theory, which is ex- 
tended to take into account the spin degrees of freedom. 
We start with a general Hamiltonian which fully con- 
siders the system's symmetries, such as translation and 
rotations in real space and spin space: 

ij 

+1 E£*M( F °)<*( F «)****']> ( 16 ) 



where 



a ijkl 



2m 



Anil ao + 2d2 



Airh a,2 — ao 



(17) 
(18) 
(19) 



m 3 

The subscripts are a — (x,y,z) and i,j,k,l = (0, ±1) 
The matrix B is a diagonal matrix whose elements are (1 



where 

Aij = h(r)8ij - Bij 



-fin 



(l^fel 2 + pkk) Sv + (Ml + p*u) 



a kl 

+ (F a ) a (F a ) kj + p kl % (24) 

^^[(fol^tftW,], (25) 

a kl 

u q (r,i) and v q (r,i) are the q-th eigenfunctions with the 
spin i and e q corresponds to the g-th eigenvalue. 

Equation (21) is rewritten in terms of these eigenfunc- 
tions as 

M r ) = £K( r '*) M ?( r 'j)/( £ 9) +U 9( r ' i ) U 9( r '-?')(/( e 9) + 1 )]- 

(26) 

We are treating a cylindrical system r = (r,9,z) and 
assume that each Boson does not have the spatial angular 
momentum. The relevant physical quantities are written 
in the following form up to the phase factor: 



</>j( r ) = 4> 3 {r), 
u q {v,j)=u q {r,j)e^ e ^, 



(27) 
(28) 
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v q (r,j) = v q (r,j)<*» e J*", (29) 
Pij( T ) = Pij(r)- (30) 

Thus the system is specified by the quantum numbers; 
the angular momentum qg and the momentum q z along 
the z direction. The density of atoms in the spin state 
i is \4>i\ 2 + pa- Therefore, the total particle number and 
the polarization corresponding to eqs. (6) and (7) are 
given by 

N = j ^(|&| 2 + Pu)dT, (31) 

i 

M = j {{\<t>i\ 2 + pxx) - {\4>-x\ 2 + p-x-i)}dv. (32) 

We have done the self-consistent calculations to solve the 
above coupled equations. The actual numerical compu- 
tations are performed under the conditions: The mass 
and the scattering length of atoms are m = 3.81 x 
10 _26 kg and ao = 2.75nm. These are equal to those 
of Na atom. The other scattering length 02 is defined so 
that g s becomes ±0.1. The area density per unit length 
along the z axis is N/L = 2 x 10 4 (/im)~ 1 . The system 
size of z direction is L = 10p.m. The trapping frequency 
is u}/(2ir) = 200Hz. The energy is scaled by the trap 
unit hio. The resulting transition temperature To for the 
non-interacting Boson system is To = 0-757pK . 

The excitation levels whose energy exceeds 50?ia; are 
calculated by ignoring the interaction energy. Most of 
the neglected energy comes form the diagonal part of 
the g n term in eqs. (24) and (25). That is approximately 
estimated as g n x (density) ~ 7Tiu>. This is sufficiently 
smaller than the cutoff bQTiuj, meaning that the energy 
levels above the cutoff are not affected by the interactions 
and the above cutoff procedure is reasonable. 

§4. Interacting Bose systems 

By taking into account the presence of the non- 
condensates, we calculate a self-consistent set of the 
physical quantities under given parameters such as g n , 
uj, T or the polarization M/N . We assume that the ficti- 
tious magnetic field B is applied to the system, which is 
treated as a Lagrange multiplier and fix the polarization 
value. Both ferromagnetic (g s < 0)) and antiferromag- 
netic (g s > 0)) cases are comparatively studied in the 
following. The results will be compared with those in 
the previous non-interacting case. We see how the sys- 
tem is affected by the interaction, in particular its spin 
dependent channel in the light of the actual experimental 
situation for 23 Na atom gases with F = 1. 

Before going into the detailed separate discussions on 
the g s > and g s < cases, we first show in Fig. 2 
the particle density (a) and the excitation level structure 
(b) for a typical example (g s /g n = —0.1 and M/N — 
0.586) at a finite temperature (T = 0.3pK). As seen 
from Fig. 2(a) the condensate density with the spin +1 
component in this case occupies the central region where 



the confining potential is weak. The non-condensate, 
which includes all the components +1,0 and -1 evenly, 
distribute rather uniformly and are pushed outwardly. 
The excitation level structure is shown in Fig. 2(b) where 
the levels are plotted relative to the ground state energy 
level situated at e = 0. It is seen that the lowest levels 
come from the spin +1 component and the other levels 
with the and -1 components are pushed upwards. 
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Fig. 2. (a) Upper panel: Density profiles of the total number (A), 
the condensate with the spin +1 component (B) and the non- 
condensate (C). Lower panel: Enlarged density profiles of the 
non-condensate. The total (solid line), the mp = +1 component 
(long dashed line), (short dashed line) and -1 (dotted line), (b) 
Excitation level scheme in the lower energy region as a function 
of the angular momentum qg at q z =0 where gs/gn = —0.1, T = 
0.3fiK and the polarization M/N = 0.586. The filled dot is for 
the levels belonging to mp = 1 and the empty box is for mp = 0. 



4.1 Ferromagnetic interaction 

Let us start with the ferromagnetic case (g s < 0). 
The condensate tends to have only one spin component 
whose spin direction is antiparallel to the magnetic field. 
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Fig. 3(a) is a stereographic view of the phase diagram 
where the relative number of the condensate is displayed. 
The curve with filled dots indicates that the particle 
number of the condensate becomes 10% of the total parti- 
cle number, signalling that the transition line on the floor 
is approached from low T . The dotted curve on the floor 
is the non- interacting phase diagram shown in Fig. 1. ft 
is seen that the T c \ line is lowered due to the interaction 
effect. A drastic change from the non-interacting case oc- 
curs at smaller N c /N and lower T region where no finite 
condensation solution is reached. Therefore, upon low- 
ering T the system first Bose condenses and at a further 
lower T it reenters the normal state instead of the second 
transition T C 2 as in the ideal Bose system. In the actual 
situations the condensation with mj? = +1 may become 
phase-separated spatially into two components mp=+l 
and -1, keeping the total polarization fixed. Since the 
system is ferromagnetic, the other spin components tend 
to be avoided. This causes the empty region at the left 
hand side of the phase diagram. 

Figure 3(b) shows the level structure of a few lowest 
levels. While the levels of wf = +1 denoted by the 
filled dots stay at the same energies, other levels (empty 
symbols) belonging to mj? =0 and -1 go down as the 
polarization M/N decreases and finally reach the e = 
line where the ground state level situates, indicating 
that the condensation with mp = +1 does not exist any 
further. This lowering tendency itself agrees with the 
discussion for non-interacting system in Sec. 2, however 
the latter does not exhibit the instability. This is caused 
solely due to the ferromagnetic interaction effect. 

4-2 Antiferromagnetic interaction 

Let us consider the antiferromagnetic case where g s is 
positive. The same set of the parameters as in the pre- 
vious Subsection except that the sign of g s is reversed, 
keeping its magnitude same, gives a quite different stere- 
ographic view of the phase diagram shown in Fig. 4. The 
phase diagram in the M/N vs T plane looks similar to 
that in the ideal Bose case in Fig. 1. There exist al- 
ways the double transitions at T cl and T c2 in the present 
case. The higher transition T c \ indicates that the spin 
component with mp = +1 appears while at the second 
transition T c2 the mp = — 1 component appears in a 
continuous manner. Note that the mp = condensate 
component never shows up. This is a different point from 
the ideal Bose system where all three components appear 
at lower T. However, the overall behavior is essentially 
the same as in the ideal Bose case. We can see slight 
depressions of T c \ and T c2 compared to the ideal case. 
This is one of the interaction effects. The magnitude of 
the T c depression depends strongly on the density chan- 
nel interaction <?„, but only weakly on the spin channel 
interaction g s as shown shortly. 

Figure 4(b) shows the excitation level structure at 
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Fig. 3. (a) Stereographic view of the phase diagram for g s /gn = 
—0.1. The axes show the relative magnetization, temperature 
and the particle number of the condensate normalized by the 
total particle number. The dotted line on the floor is T c \ and 
T c 2 lines of the non-interacting system, (b) Excitation levels as 
a function of M/N at T = 0.3fiK. The solid line means the 
angular momentum qg = and the dashed line means qg = 1. 
The filled dot means levels belonging to mp = 1. The empty 
box means mp = and the empty triangle means mp = — 1. 



the lower energy region. Here the levels belonging to 
mp = —1 and mp = never reach the e = line as 
M/N decreases. There is a zigzag of the levels at around 
M/N = 0.5. This signals the second transition point at 
T c2 , corresponding to a = in §2. 

In order to examine the interaction effects quantita- 
tively, we also calculate another set of the parameters 
for the system having a weaker g n constant, compared 
with the above result shown in Fig. 4 while keeping the 
ratio g s / g n = 0.1 fixed. The two transition temperatures 
differ from those in the ideal Bose case discussed in §2 
and also slightly from that in Fig. 4. The resulting phase 
diagram is in between them, namely, the two transition 
temperatures become similar to that of the ideal Bose 
case. 

§5. Conclusion 

In this paper, we have examined the phase transition 
in the spinor BEC with the hypcrfine state F = 1 . We fo- 
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Fig. 4. (a) Stcrcographic view of the phase diagram for g s /g n = 
+0.1. The axes show the relative magnetization, temperature 
and the particle number of the condensate normalized by the 
total particle number. The dotted line on the floor is T c i and 
T c 2 lines of the non-interacting system. The dashed line means 
the mp = +1 fraction of the condensate. At smaller m, the —1 
fraction of condensate appears (between the dashed line and the 
upper solid line), (b) Excitation levels as a function of M/N 
at T = 0.3fiK. The solid line means qg = and the dashed 
line means the angular momentum qg = 1. The filled dot means 
levels belonging to mp = 1. The empty box means mp = and 
the empty triangle means mp = — 1. 



tially. 

The interaction g n of the density channel acts simply 
to suppress both transition temperatures T c \ and T C 2 by 
similar amount. Therefore it is found that g n modifies 
the system quantitatively while g s drastically alters in a 
qualitative manner. 

We also study the spatial distribution of each spinor 
BEC component and the excitation spectral structure. 
It is our hope that by using optical trapping method 
these predictions must be checked experimentally since 
we have the antiferromagnetic spinor 23 Na BEC and 
87 Rb is expected to be ferromagnetic. 
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cus on the interaction effects: There are two kinds of the 
interaction channels in the present spinor BEC, namely, 
the spin channel specified by g s and the density channel 
specified by g n . The ideal Bose system always exhibits 
the double transitions from the normal state to the one- 
component BEC at T c i and then to the three component 
BEC at T C 2 upon lowering T. Depending on the nature 
of the interaction, or the sign of g s , only the antiferro- 
magnetic spinor BEC (g s > 0) shows the double tran- 
sitions. In this case the nif = -1 component becomes 
non- vanishing below the lower T c2 where the mp = +1 
component appeared at the higher transition T cl persists 
in the ground state. 

In the ferromagnetic case (g s < 0), on the other hand, 
the spinor BEC always shows the instability at a lower 
temperature. This lower instability transition may actu- 
ally be realized as the phase-separation. The dominant 
rriF = 1 component and other components segregate spa- 



